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DAMPING IN ROLL OF CRUCIFORM AND SOME RELATED 
DELTA WINGS AT SUPERSONIC SPEEDS 
By Herbert S. Ribner 


SUMMARY 

The damping in roll of cruciform delta wings in supersonic flow 
has been evaluated by means of small-disturbance (linearized) wing 
theory; both subsonic and supersonic component stream velocities 
normal to the leading edges have been considered. In addition, some 
known two-dimensional results for rotating multibladed laminae have 
been applied to obtain the loading when the number of panels is 
changed from four to an arbitrary number, under the restriction of 
low aspect ratio:, the damping in roll has been determined explicitly 
for three panels. Finally, the damping for an infinite number of 
panels has been evaluated without restriction as to aspect ratio or 
Mach number. 


INTRODUCTION 

The linearized theory of planar wings in supersonic flow has 
been highly developed. The theory for nonplanar wings such as the 
'cruciform seems to have been limited, until very recently, to the 
analysis (reference 1) of the cruciform delta wing with angle of 
attack and. sideslip. The rolling behavior of a cruciform delta 
wing is likewise of considerable importance in flight mechanics. 

The present paper is concerned primarily with the evaluation of the. 
damping in roll. 

The cruciform delta wing is illustrated in figure 1(a). This 
type of wing develops a large side force when sideslipped, thus per- 
mitting sharp turns without banking. With or without sideslip, the 
problem of calculating the lift distribution is simple because the 
horizontal and vertical wings behave as though they were isolated. 

For rolling motion, however, the situation is more complex. 

The horizontal and vertical wings mutually interfere, so that the 
pressure distributions are not the same as if the wings were iso- 
lated. The calculation is no longer a planar problem, and the well- 
known planar techniques in linearized wing theory must be extended. 
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For the very narrow (low-aapect-ratio) cruciform, the Munk- Jones 
approximate theory (reference 2) may be applied. The problem 
reduces to the determination of the potential of a rotating cruci- 
form lamina in two-dimensional flow. This problem has already been 
solved in connection with propeller theory; in addition, the prob- 
lem of the rotating lamina has been solved (in series form) for 
any finite number of blades (reference 3). This work is discussed 
and applied in the section "Low-Aspect-Ratio Approximation for Multi- 
planar Wings." 

The case of a broad cruciform wing with supersonic leading 
edges may readily be reduced to a planar boundary-value problem. 

The principles of superposition and reflection are used. When the 
leading edges are subsonic, a partial reduction to a planar problem 
is still possible. The solution is then completed by an iteration 
procedure of rapid convergence. The simpler supersonic- edge case 
is treated in the section "Cruciform with Supersonic Leading Edges" 
and the subsonic-edge case is treated in the section "Cruciform 
with Subsonic Leading Edges." 

Two contemporary investigations covering portions of the scope 
of the present investigation have recently been noted. The slender 
cruciform has been treated by Gaynor J. Adams of the NACA Ames lab- 
oratory. The cruciform delta with supersonic leading edges is 
reported in reference 4. The cruciform delta with subsonic lead- 
ing edges, which occupies most of the present investigation, 
apparently has not been treated elsewhere. 

LOW-ASPECT-RATIO APPROXIMATION FOR MULTIPLANAR WINGS 
Arbitrary Number of Panels 

The multiplanar wing under consideration is illustrated in 
figure 1(b). Each panel is an identical right triangle, and the 
long .leg is common; the cross section perpendicular to this leg 
presents a spoke-like appearance. The damping in roll about the 
axis of symmetry is to be calculated. 

The restriction to low aspect ratio (triangular panels very 
slender) allows the approximate theory of reference 1 to be 
employed; the results apply for subsonic, transonic, and a range 
of supersonic speeds. According to this theory the local potential 
at any section perpendicular to the axis of symmetry is approxi- 
mately the same as if the section were a two-dimensional lamina. 

The present problem is thus reduced to the determination of the 
potential distribution on a two-dimensional lamina of N equally 
spaced spokes, rotating about the center. 
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It happens that a similar reduction of the problem of the wake 
of a multlblade propeller Is possible. As the wake pitch approaches 
infinity, that problem likewise reduces to that of the two-dimensional 
rotating lamina. In this connection, the potential distribution for 
the lamina was obtained in reference 3; there is evaluated therein 
a function K that is related to the surface -potential Jump A <f> 
by 


A0 = 2«JLP k 
v N 

in the present notation. This equation, together with the equa- 
tion for K in reference 3 (p. 682), gives the series for the 
potential Jump 


a 0 = 2 2 - 4 / n s 2 p rq-ti/NL,, 

[r (l+2/N)] 2 ^ 


(-l) n+1 (N-2) (2N-2) . , . (nN-N-2) ain Kn0 
(N+2) (2K+2) . . . (nN+2) 


( 1 ) 


where the parameter 0 is defined by 


r = s 


(“ \ *>) 


2/N 


and r is the radial distance. The series "is troublesome to 
handle analytically but is quite convenient for numerical 
computat ion . " 


Three Panels 

The insertion of N = 3 reduces equation (1) to 


[r(2/3)J 2 n=l 


( 2 ) 


with r = 3 


|cos ( 30 / 2 jj 


2/3 
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The series in equation (2) has been summed, numerically in refer- 
ence 3. The results have been used, to obtain the values of 

MA0 - , ■ 

— — against r/s tabulated in table I. 


The rolling moment developed by the three panels is given by 


L = -3pV / A0 r dr 

J 0 


The integration has been performed graphically and the result is 

L = - 0.531 pVps 4 ( 

with about 1 percent uncertainty. 

This damping compares with that (reference 5) for the flat, 
narrow delta (2 panels) as follows 


damping in roll, N = 3 
damping in roll, N = 2 


= 1.35 ± 0.014 


This ratio is about 10 percent less than the ratio 1.5 that would 
be obtained from an elementary calculation neglecting interference. 
The reduction is caused by mutual interference among the three 
nanels . 


Four Panels (Cruciform) 

The insertion of N = 4 reduces equation (l) to 


A0 = - s 2 p ^ (-1) 
n=l 


n+1 sin 4n0 
4n 2 -l 


This series may be summed, and according to reference 3 the result is 


a a 2 2 0 . , 1+tan 0 

A0 = - a p cos 20 log Q 

it e 1-tan 0 


( 6 ) 
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with r = s *Jcoa 20. The span loading is proportional to 

NA0 

ths variation of with r/s is given in table I. 

a P ' 

The rolling moment developed by the four panels is given by 


1 


L = - 4pV / A0 r dr = - 2p ^ ~ 


(7) 


Comparison with the flat, narrow delta 

damping in roll, N = 4 
damping in roll, N = 2 

If there were no mutual interference between the panels, the factor 
would be 2. The loss in damping caused by the interference is thus 
19 percent for the narrow cruciform delta. 


(two panels) gives 
16 

= — = 1.62 


( 8 ) 


The nondimensional damping- in-roll derivative C 7 for the 

P 

narrow cruciform is obtained by dividing equation (7) by pVps*0/m. 
The resulting coefficient, which is based on the area of Just two of 
the four panels, is 


C 



(9) 


Infinite Number of Panels 

Squation (l) is inconvenient for evaluating the behavior when 
the configuration has an infinite number of panels (N =<»). The 
result can be obtained, however, from simple physical considera- 
tions, which need not be limited to the case of a narrow wing. Fix 
attention on the section perpendicular to the axis in figure 2, 
and imagine the number of panels greatly increased. Qualitatively, 
the air trapped between the panels will tend to move as a unit with 
the blades; so will the 'cylindrical wake. (The centrifugal forces 
acting on the rotating air mass are second-order effects proportional 
to p^ and hence may be neglected in the present analysis.) Fur- 
thermore, the external disturbance will be like that near the edges 
of a moving cascade: this disturbance will be appreciable only 
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within radial distances comparable with the gap distance between 
panels at the periphery. In the limit , then, as the number of 
panels becomes infinite, the external disturbance vanishes, and 
the interior fluid (that within the conical envelope of the panels 
together with the cylindrical wake) rotates as a solid. 


This result may be proved as follows : Assume that which 

is equal to the circulation F, varies as r^ in the radial direc- 
tion along the panel and is independent of x. This implies that 
the loading is concentrated along the leading edge and that dF/dr 
varies as r. The total circulation threading an annulus of radial 

thickness dr, coaxial with the multiplanar wing, is ^ 7 - dr. The 

mean vorticity in this annulus is the circulation divided by the 


area 2nr dr. This vorticity is , which is independent of r. 

2 rtr 

As the number of panels N approaches infinity the vorticity 
approaches a uniform distribution within the annulus. Because there 
is no change from one annulus to another (independence of r), the 
vorticity is uniform over the entire circular cross section normal 
to the axis. This property holds for all such cross sections from 
the apex of the multiplanar wing back into the cylindrical wake. 


Such a uniform distribution of vorticity implies that the fluid 
bounded by the conical envelope of the wing and the cylindrical sur- 
face of the wake rotates as solid. The outer fluid is completely 
undisturbed, because a surface sheath of vorticity of reversed sign 
achieves exact cancellation of the disturbance farther out. This 
surface sheath results from the discontinuous drop of F to zero 


at the edge of each panel 


- - -). 


ylr 


The rotation of the interior 


fluid as a solid clearly complies with the boundary. conditions at 
the panels. Moreover, the Prandtl-Glauert partial differential 
equation is satisfied between the panels and in the outer space. 
Thus the assumed parabolic distribution of A0 = T (which may be 
interpreted as a distribution either of doublets or vortices) satis- 
fies the boundary-value problem and must be the correct solution. 


The damping in roll may now be evaluated either by considera- 
tions of angular momentum or of individual panel loading. The 
momentum approach is given because of its greater simplicity; the 
panel-leading approach gives identical results. The torque or 
rolling moment is equal to the time rate of increase of the angular 
momentum of the wake. This may be written 


L = - IpV 
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cylindri- 
for I 

( 10 ) 

Comparison of this damping moment with that for a flat, narrow 
delta (N - 2) gives the proportion 

damping in roll , N = op _ 4 
damping in roll, N = 2 1 

In other words, the addition of an infinite number of similar panels 
radially disposed to a flat narrow delta wing increases the total 
damping in roll by just a factor of 4. 

The foregoing analysis for an infinite number of panels does 
not involve the assumption that the panels^ are slender nor any 
restriction on the Mach number. Thus the damping- in-roll equa- 
tion (10), unlike the corresponding equations for three and four 
panels, is valid for all aspect ratios and any speed, subsonic or 
supersonic, so long as the disturbances are small. 


where I is the moment of inertia per unit length of the 
cal rigidly rotating wake. Insertion of the value ps 4 
yields 

L = - | pVps 4 


CRUCIFORM WITH SUPERSONIC LEADING EDGES 

Assumptions. - The assumptions of an inviscid fluid and small 
disturbances (isentropic flow) have been made throughout this paper. 
Thus the disturbance velocities u,v,w possess a potential 0 
that satisfies the linearized equation of motion 

(l-M 2 )^ + 0yy + <f> Z Z » 0 ( 12 ) 

known as the Prandtl-Glauert equation. To the first order, the 
relation between pressure and velocity disturbances is 

P = - pVu (13) 

where 


a . t x 
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Analysis . - The boundary conditions for the cruciform wing 
rolling clockwise with angular velocity p are shown in figure 2, 

The conditions are 

w - - py 

on. the horizontal panels, and ■ . . 

v = pz 

on the vertical panels; in addition, u, v, and w must vanish out- 
side the envelope of the Mach cones from the apex and the leading 
edges. 

t 

The flow in each quadrant is the same, except for orientation, 
because of the fourfold symmetry. It will thus suffice to limit 
attention to the first quadrant, defined by AOB. Because of the 
linearity of equation (12), the principle of superposition may be 
used. Accordingly, the flow specified by the boundary conditions 
on AOB (redrawn in fig. 3) may be obtained by the linear super- 
position of the flows defined as 3(b) and 3(c). 

\ 

The general character of the flow specified by the boundary 
conditions labeled Flow 3(b) may be anticipated. The value of u 
will differ from zero on both OA and OB. Call the value on OA 
(upper surface) u A , and the value on OB (right-hand surface) Ug. 
New flow 3(c) is just a mirror image of flew 3(b) rotated through 
90° with sign change. This flow is characterized by a u-velocity -u A 
along OB and a u-velocity -Ug along OA. Superposition yields 
the resultant values: 

£u = (Iu) A = u A - u B (14) 

along the upper surface of OA, and 

Lu = (Zu) B = - (u^) (15) 

along the right-hand surface of OB. 

The problem of calculating the pressure distribution on the 
rolling cruciform wing (fig. 2) has now been reduced (equa- 
tion (13)) to the problem of calculating ua and ub for flow 3(b). 
This problem is easily solved by the method of images. The method 
is applicable because the vertical panel OB completely isolates 
disturbances from the left and right sides. Thus consider fig- 
ure 4(b). The right-hand side is identical with the figure speci- 
fying flow 3(b), and the left-hand side is a mirror image, the 
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combination being symmetrical about OB. The symmetry automatically 
satisfies the condition v = 0 along OB. By virtue of the sym- 
metric linear variation of angle of attack, the delta wing COA is 
called a wing of linear twist or linear washin. 

Solution for u^. - The calculation will employ figure 4(a). 

Regions I and II, separated by the Mach cone from the apex, must 
be distinguished. Region I is uninfluenced by the remainder of 
the wing, and the pressure distribution and u-velocity there are 
the same as for an Isolated rolling delta wing. The value for u 
therein may be taken from reference 6, table II, or reference 7. ; - • 

In the present notation it is 

u . gft fehU - region I (16) 

■ e 2 (m 2 -l) 3/2 

where O = Py/x 


The solution for the u-velocity in region II may be effected 
•by standard methods (for example, reference .8). The result is 

Pq2z - 3 -7 3 [- 2^/(m 2 -l)(l-<£) + (1-mo) cos -1 + 


U A,II = 


rtP 2 (m 2 -l) 


(1+mo) cos - ^ l+ mCT l 

m+a J 


(17) 


Solution for Ug. - Equation 3 of reference 8 gives the sur- 
face u-velocity in terms of a surface integral of dw/dx and a 
line integral of the value of w along the leading edge. The for- 
mula is easily generalized to values of u above the surface. For 
the present case, only the line integral survives and 

«B = V x *°> z ) 




w dy x 

\2 l 2 2 „2 2 

yix-xj) -p y-L -p z 
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On the wing w = -p kil and along the leading edge 8 |y^| = mx^. 

Because of the symmetry, the integral can be doubled and limited to 
the right-hand leading edge. The substitutions yield 


"B 


_ 2pm 



*i di i 


n P 0 ‘^/(x-x 1 ) 2 - nrx-^-p 


2 2 
z 


where M.C. denotes the value of x-^ on the Mach cone, where the 
square root vanishes. Finally 


^ = 


2>f f X 3/2 lyt- 2 - 1 ’^- <?> - tan^^-lXl-a 2 ) (18) 


«p (m -1) 


where o has been written for (3z/x. If, now, the earlier signifi- 
cance Py/x is restored to O, equation (18) may be applied along OA 
(fig. 3(a)) as desired. When the leading edges are sonic, m = 1, 
and by application of a limiting process, equation (18) is seen to 
reduce to 


-B = 5EL (l-o 2 ) 
3rtp 


3/2 


(19) 


Resultant u- velocity and pre ssur e distribution . - In region I 
Ug vanishes; therefore, according to equations (l4j and (16), 


(Zu)i = u A ^ 


pm 2 x (mo-l) 
3/2 

eW-i) 


region I 


( 20 ) 


For region II, substitution of equations (17) and (18) in equa- 
tion (14) gives 
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(lu) 


ii " — = gq2 ^ "5 72 (r 2 tan ‘ 1 ^ m 2 * 1 H 1 -° 2 ) + 

jtp 2 (m 2 -l) 


(1-niO) cos” 1 1 ‘*° 1<J + (l-Hno) cos” 1 1 ~ K ^° 
m- o m+O _ 


region II 
( 21 ) 


Equations (20) and (21) give the resultant u-velocity on the upper 
surface of OA. By equation (15) the value on the right side of 
OB Is the same thing with opposite sign (with a designating 
pz/x instead of Py/x as for OA). 


The surface pressure is related to the velocity by equation (13), 
Thus pressure is obtained by multiplying the values in equations (20) 
and (21) by -pV. 


Damping in roll, - The pressure distributions just determined 
for each of the two inner faces of the first quadrant (fig. 2) 
clearly contribute identical additive rolling moments. The total 
rolling moment for the cruciform is therefore eight times that 
contributed by the upper face of OA. 


The local pressure on OA is proportional to lu and by equa- 
tions (20) and (21) Zu is of the form x f(o). For such a dis- 
tribution calculation shows that the rolling moment contributed by 

■*m 


this pressure is proportional to 


of(o) do. The damping- in- 


roll coefficient for the complete cruciform is proportional to this 
integral and is given by 



where the factor of eight mentioned previously has been taken into 
account. Upon substituting equations (20) and (21) and carrying 
out the integration, there results 



‘.'j ’ , i j . r 






5.V/V 

ht 


m- 


0: 


W- 


m i 
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3rtp(m 2 -!) 2 


2m 2 -5 + an 4 j:5 m 2 46 , gee' 1 m 


^m 2 -l 


Equation (22) gives the damping- in-roll coefficient for a 
cruciform delta wing with supersonic leading edges. For the case 
of sonic leading edges (m = 1) a limiting process yields 


(This same result can he obtained with considerably less labor by 
limiting the foregoing analysis at the outset to the sonic-edge 
case . ) 


CRUCIFORM WITH SUBSONIC LEADING EDGES 


The cross section of a single quadrant in a plane x = constant 
is shown in figure 5. The trace of the Mach cone from the apex is 
also shown. The wing panel traces are shown solid, and dotted-line 
extensions out to the Mach cone have been added. On the two wing 
panels the boundary conditions for rolling motion (angular velocity 
are w = -py and v = pz, respectively. The dotted-line exten- 
sions carry no lift, and this condition can be specified by the 
relation Au = 0, Au being the jump in u across the dotted line. 
Finally the entire disturbance u,v,w must vanish on the Mach 
cone. 


It is possible to satisfy these boundary conditions by means 
of a suitable distribution of doublets on the panels, the same on 
each panel. When the fourfold symmetry is taken into account, the 
u-velocity induced by these doublets at the dotted lines is seen 
to be zero. Thus the conditions Au = 0 become u = 0. With this 
change then, the boundary conditions for adjacent quadrants are no 
longer interdependent and the flow in the typical quadrant (fig. 5) 
is completely determined by the conditions along its own boundary. 


Figure 5 is redrawn' with u = 0 in place of Au = 0 and des- 
ignated flow 6(a) in figure 6. Starting with this figure, the 
problem is further simplified by the principle of superposition. 

The flows specified as 6(b) and 6(c) will superpose to yield the 
flow specified as 6(a). Flow 6(c) can be seen to be a mirror image 
of flow 6(b) rotated through 90° with sign change. The problem for 
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the entire cruciform is thus reduced to the solution for the flow 
specified by the boundary conditions labeled flow 6(b). 


Because of the difficulty of obtaining an exact solution, an 
iteration procedure was decided upon. The procedure is schemati- 
cally represented by figure 7. The unknown u-velocity in the 
specifications for flows 7(a), 7(b), 7(c), . . .is to be solved 
for successively. The flows may then be superposed to converge to 
a solution for flow 6(b). Flow 7(a) clearly violates the require- 
ment u = 0 in the vertical dotted region. Flow 7(b) attempts to 
rectify the situation by cancelling the velocity u = ug obtained 
in that region. In so doing, however, flew 7(b) leads to a much 
smaller induced u-velocity in the horizontal dotted region. This 
velocity u = u A ^ is cancelled by flow 7(c), but at the cost of 

again slightly violating the condition u = 6 on the vertical 
dotted region, and so on. 


If the dotted regions are relatively small (parameter m not 
too small compared with unity) the successive u- velocities in 
flows 7(a), 7(b), 7(c), . . . , will diminish rapidly, providing 
rapid convergence. It will be shown that the iteration process may 
be terminated after the evaluation of flow 7(b) with small loss in 


accuracy if m > 0.5; the velocities to be calculated are then 


"B' u B,b> 


Solution for u A . - The vertical panel bounding flow 7(a) 

(fig. 7) completely isolates disturbances from the left and right 
sides, and the method of images is applicable, as for the supersonic- 
edge case. The same reasoning as for that case leads to figure 8, 
which is again a delta wing with linear twist, but this time with 
subsonic leading edges. The value of u on the right-hand panel 


of this wing (upper surface) is the desired value of 


The solu- 


tion presents a problem in its own right, and the details are given 
in appendix B. (The solution was also obtained in reference 9 by a 
different method.) The result is given in equations (B2) and (B17); 


u A = 2 ? x a(m) + ^ b(m) 

4FJ 


where 
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a(m) , 3m 5 (E-m 2 K) 

(l-2m 2 )E+m 2 K 


b(m) = 5m5 ^- E ). 


(l-2m 2 )E+m 2 K 


and K and E are complete elliptic integrals of modulus k 


= 1-m 2 . 


Solution for 


U B* 


- According to reference 10 the analytic 


continuation of a conical surface distribution of u, u( 0) into 


space is obtained (subject to certain restrictions) by replacing O 
by £ and taking the real part. (See discussions in appendix B.) 
In the present case u^ is not conical, but du^/c>x is conical, 

and the procedure will apply to the latter. The equation for 
du^/dx with £ inserted is given in equation (B7). Along the 

z-axis (that is, along OB) £ assumes the value iu/^/l-u 2 where 


V = Pz/x. (See equation (B3).) With this value equation (B7) 
becomes 



2P 


a m 2 +ey 2 


B 3rtp 


[m 2 +(l-m 2 )u^] 


z/2 




where 


e = (2-m 2 )a + m 2 b 

The integration with respect to x then yields 


«b - f,<*> 

3rtp 


where 


f B ( u ) 5 


5 e + -^St) + 



and 












k = y/l-m 


(This expression reduces to the simple equation (19) when the lead- 
ing edges are sonic, upon applying a limiting process as m — fl.) 
The variation of Ug with u is shown in figure 9 for m = 0.866. 

Solution for Ug g. - Refer again to flow 7(b) (fig. 7). The 

condition w = 0 along the horizontal axis permits the use of 
reflection in this axis. Thus Ug g is the value of u on the 

right-hand panel (lower surface) of the wing shown in figure 10. 

The exact solution of this wing presents an almost insuperable prob 
lem, and even a close approximation (for example, by the method of 
reference 11) is unduly laborious. The average magnitude of Ug g 

can be anticipated, however, to be small enough compared with Ug 
and u A so that some inaccuracy may be tolerated. Accordingly, in 
place of the exact Ug g a relatively easily calculated flow will 

be used that fails to cancel ug exactly for u > n, but does can- 
cel it on the average. 

This approximate distribution Ug g may be specified with the 

aid of figure 10 with z replaced by y and upper and lower sur- 
faces interchanged for convenience. This flow consists of the 
superposition of four flows, (l), (2), (3), and (4), determined by 
the following respective boundary conditions at the wing surface: 


(1) -1 < V <1 

w 

= -Kx 

(2) -1* U s>l 

w 

= +*P|y| 

(3) -m < u < m 

w 

= KX 

m < M < 1 

u 

= 0 

(4) -m< u< m 

w 

= -Kp|y| 


m < luj <, 1 u = 0 


It is clear that the superposition of these flows satisfies the 
required boundary condition for Ug g that w = 0 in the 
region -m < U <_m. In addition, t&ere will be a nonvanishing 
quasi-conical distribution of u in the regions m < ^I<J1 from 
flows (1) and (2), although flows (3) and (4) contribute nothing 
there . 



I < *4' 


<^*i;«^V..iJ V>»-' f ;:y .^«?3>v" ^v? *" 
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Elows (1) and (3) may be recognized as corresponding to the 
pitching delta wing with sonic and subsonic leading edges, respec- 
tively. Flows (2) and (4) correspond to the delta wing with linear 
twist, again with sonic and subsonic leading- edges, respectively. 
The solutions for (l) and (3) are well-known (reference 12) and 
those for (2) and (4) have already been obtained herein. The solu- 
tions for the u-velocities are 


4 toe 2- 


^7 


-2KX l+l? 
3«P . / 


V 7 F 

2kx «i (m)+ ^ b l (n) 


/m 2 V 


U 4 = 3rtP 


2kx a(m)+u b(m) 


/ 2 2 
fm - u 


where a^m), b-^m), a(m) , and b(m) are defined in appendix A. 
Note that u^ and u^ are to be taken equal to zero for v > m. 

The sum u i + ^ + ^ + u 4 is " fche ten ‘ ba ' fcl - V9 approximation to 
Ug This sum is 


^b " 


l_ u 2 _ A(m)-t- iTB(m) 




-m < u < m 


(26a) 


m < lu si 


(26b) 


where 


A(m) s a 1 (m) - a(m) 


B(m) = b^m) - b(m) 

and the parameter * is yet to be determined. 
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In the region m < |vi <_ 1, Ug g is given by equation (26b). 

The variation with U ia shown graphically in figure 11 for a 
selected value of Cn the same figure ia shown a portion of 

the variation of Ug with u for m = 0.866. The value of 

was chosen to make ug g = -0.8 Ug at v = m+. This value is, 
by virtue of equations (25) and (26b) , 


0.8 f £ (m) p -0.8 f B (jn) p 
/ f 3 = 3k 8 

3Vl-m 


(26c) 


where fg(m) is fg(u) (equation (25)) with u = n. This select iai 
of k doe 8 not effect an exact average cancellation of Ug by 
Ug g in the range m< |bl<,l; instead, the ordinates near V = m 

are given more weight as suggested by equation (17a) of reference 11. 
It is thought that the average of Ug g calculated in this way 

over the range 0 £. v <m is within ,±25 percent of the exact value. 


Solution for u A g. - The distribution u A g along the y-axis 

(fig. 7(b)) is obtained by analytic continuation of the distribu- 
tion Ug g, which is known along the z-axis. The procedure is pat- 
terned after that discussed in the solution for ug. The result 


is obtained without difficulty as 


u 


A.b 



k) - 2E(0,k) + 


E(0,k) - nrF(0,k) - 


A 

's/m 


2 + k 2 o2 h. 


where 


(27) 


eg = (2-m 2 ) A + m 2 B 


Resultant u-velcclty and pressure distribution. - The unknown 
surface u-velocity distribution in the flow defined by figure 6(b) 
has now been determined to the desired accuracy. The value on the 
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upper surface of QA is u A (o) (equation (24)) plus u A b (o) (equa- 
tion (27)); the value on the right-hand side of OB is Ug(u) 
(equation (25)) plus Ug g(u) (equations (26a) to (26c)). To the 

flow of figure 6(h) must be added the mirror- image flew with 
reversed sign defined by figure 6(c). This image flow is char- 
acterized by a u-velocity -u A (u) -u A ^(u) along OB and a 

u-velocity -Ug( d) ' -ug g(a) along OA. Superposition yields the 
resultant values: 

Eu = (Z U ) A = u A (o) - ug(o) - Ug^Cd) + u A>b (o) (28) 

along the upper surface of OA, and 

Eu =.(Zu) B = - [u A (w) - Ug(u) - Ug^(u) + u Ajb (y)j (29) 

f 

along the right-hand side of OB. 

The surface pressure distribution (deviation from stream pres- 
sure) is obtained by multiplying the values of Zu by -pV, 
according to the linearized Bernoulli equation (13). ' 


Damping in Roll 


The pressure distributions just determined for each of the two 
inner faces of a quadrant of the cruciform contribute identical 
additive rolling moments; this is a consequence of the anti- 
symmetry of equations (28) and (29). The total rolling moment for 
the cruciform is therefore eight times that contributed by the 
upper surface of OA. 


The local pressure on OA is proportional to lu, and the 
equations show that each of the components of £u is of the 
form xf^(o) . For each such component the contribution to the 

damping- in-roll coefficient for the entire cruciform was earlier 
pointed out to be 


A i C 2 



of i (o) do 
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Thus the total value of Cj is 

V 


C Z p = A A C Z p + A B C 2 p + A B,b C I p + A A,b C l 


(30a) 


the increments corresponding respectively to u A , -Ug, -Ug B an ^- 
u A i^. Upon carrying out the integrations there is obtained 


4 

9«pm3 


|j3ma + 2m 3 hj 


(30b) 


±- m\(m) +‘i (l 


9rtPm 


h) v ■ y + 2 ^ 


-1 ,.2 



k 3 


_ 4 

0.8 fg(m) 

fraA + — 

9jtPm 3 _ 

k 

V. 3 

4 

0.8f B (m) 

d 

r 3 

l+3k-4k° 


(30c) 


(30d) 


‘ A > b 


:'\/l+k 2 


n- 


. k — A (a -!f|+ sln lk - 8ln ^ L 

2k 1 * 3 L 


(3-2k 2 ) e l 
2k 2 

(30e) 


vhere a, b, e, A, B, and e-^ are functions of m defined in 
appendix A; the function f B (m) has been evaluated from f B (h) in 
equation (25) 
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f B (m) - | (- m 2 F(^,k) + 


(- m 2 F^k) + E(^,k) - 


(30f) 


where 

^ = sin ' 1 k 

and f A - 5 ( 111 ) is fg(m) with A in place of a and B in place 
of b. 

The variation of -C7 with m according to equations (30a) 

P 

to (30f) has been plotted in figure 12. On the same figure the 
successive contributions of the iteration process, A A c l > a B c l > 

P P 

Ar> V C, , and A. V C, are shown. It is seen that the sum A„ ,C, 

B > B V *p *p 

plus A A jjCj reaches only 4.7 percent of Cj as m is reduced 
r P P 

to 0.5 and thus justifies the roughness of the approximation for 
Ug -5 and u A - 5 . The estimated +25-percent uncertainty in the aver- 
age value of Ug -5 and u A 5 (which is surely liberal) accord- 
ingly contributes only ±1. 2-percent uncertainty to C 7 at m * 0.5, 

*P 

A plot of % 5 (hot shown) shows a general similarity to 
that for ug (fig. 9). From this similarity the contribution to 

C. of the next iteration, flow 7(c) of figure 7, with velocity 
P 

components u A ^ c and ug^ c , can be roughly estimated. The pro- 
portion contributed is according to this estimation: 


A C7 +A Ci 
A.c *p B.c 1 


-u (0)\/A C 7 +A C 7 

aV Y B,b ! p A.t 1 

V o ) A °i„ 


This ratio increases with decreasing m; for m = 0.5 the value 
is (0.025)(-0.047) s -0.0012, or a little over -0.1 percent. Thus 
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above m = 0.5 the total contribution of flow 7(c) is small com- 
pared with the uncertainty (±1.2 percent of C 7 at m = 0.5) in 

L V 

the contribution of flow 7(b). For this reason flow 7(c) and fur- 


ther iterations have been neglected in the computation of Cj . 


Examination of figure 12 in the light of the foregoing con- 


siderations shows that the estimated uncertainty in C 7 increases 

•ip 


progressively from zero at m = 1 as m is decreased. The 
graphs are terminated at m = 0.5 where the error reaches an esti- 
mated ±1.2 percent +0.1 percent. 


DISCUSSION AND CONCLUDING REMARKS 


In the preceding three sections, the damping in roll for the 
cruciform delta wing has been evaluated for three speed regimes 
characterized by three ranges of the edge-slope Mach number param- 
eter m. Equation (9) (rigorous fcr m — >0) is considered to be a 
good approximation for 0 < m <. 0.25. Equation (30) was shown to 
be a good approximation (progressively better with increasing m) 
for 0.5 < m <, 1. Finally, equation (22) is exact for 1 <_ m < 
(Both equations (22) and (30) reduce to equation (23) for m = 1.) 


The variation of -3Cj with m from these three equations 


is plotted in figure 13. It is seen that equation (30) appears to 
be asymptotic to the slender-wing- theory result (equation (9)) for 
decreasing m. In fact, a satisfactory interpolation (shown dotted) 
in the neglected range 0.25 < m <.0.5 is indistinguishable from 
equation (9). 


The damping in roll of a multiplanar delta wing of N panels 
is always less than N/2 times the damping of the plane delta 
wing of two panels. The reduction is caused by interference 
between the panels in different planes. The present results show 
that the loss fcr the cruciform delta varies from 19 percent as 
m— +0 to 7 percent at m = 1 and decreases to zero as m — »<*>. 

For m — >0, the interference loss increases progressively with the 
number of panels as follows: three panels, 10 percent; four panels, 

19 percent; infinite number of panels, 100 percent. . The relative 
total damping in roll for m — »0 is: two panels, 100 percent 

(standard of comparison);, three panels, 135 percent; four panels, 

162 percent; infinite number of panels, 400 percent. Finally, the 
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damping in roll for an infinite number of panels, ‘regardless of 
aspect ratio, is independent of Mach number both below and above 
the speed of sound. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, August 24, 1950. 
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APPENDIX A 


SYMBOLS 


The following symbols, are used in this report: 


A = A(m) = 


a = a(m) = 


a^(m) - a(m) 

5m 5 (E(k)-m 2 K(k)) 
( l-2m 2 )E( k) +m 2 K ( k ) 


= a 1 (m) = 


3nm^k^ 


(l-2m 2 )E(k)+m 2 K(k) 


B = B(mj = 


b^(m) - b(m) 


b = b(m) = 


5m 5 (K(k)-E(k)) 

(l-2m 2 )E(k)+m 2 K(k) 


b x = b-^m) = 


- 3*^/2 


(l-2m 2 )E(k)+m 2 K(k) 


damping- in-roll derivative = 
rolling moment 


— pV^ ‘area two panels -22 


on t, an tj dn t Jacobian elliptic functions 


E = S(k) 


complete elliptic integral of second kind, 
modulus k 


E(0,k) 


incomplete elliptic integral of second kind, 
modulus k, amplitude 0 


e = e(m) = 


(2-ra 2 ) a(m) + m 2 b(m) 


(2-m^) A(m) + m^B(m) 



mm m 
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incomplete elliptic integral of first kind, 
modulus k, amplitude 0 


function of a (different definitions in dif- 
ferent places) 


f-D(m) with A in place of a and B in 


place of b 


[f B (*)U given by equation (30f) 


defined in equation (25) 


moment of inertia 


function related to A0 defined in equation (l) 


complete elliptic integral of first kind, 
modulus k 


modulus of elliptic integrals 


rolling moment, positive in sense of a right- 
hand screw proceeding upstream 


Mach number = 


flight speed 


sonic speed 


value 0 at right edge of delta wing 


/ tan c j 
\tan nj 


number of panels (half-deltas) in multiplanar 
wing (N = 4 for cruciform) 


local pressure minus stream pressure 


angular rolling velocity, positive in sense of 
a right-hand screw proceeding upstream 


angular pitching velocity, positive in sense of 
a right-hand screw proceeding outward along 
y-axis 


real part of 


HVOP 
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1 

radial distance from x-axis 

a 3 

semispan 

9 t 

argument of elliptic functions (=F(0,k)) 

i 3 u(0 

•£« CO 9 

function of complex parameter C 

m u > w 

disturbance velocity components along x-, y-, 
z-axes, respectively 

I v 

free-stream velocity (flight speed) 

m x > z 

Cartesian coordinates: x-axis parallel to free- 

stream direction; y-axis horizontal and toward 
right, looking upstream; z-axis vertically 
upward 

p p 

Vm 2 -i 

il ' ^ 

circulation 

1 r(N) 

gamma function (equation 481, reference 14) 

j j ■ ' A 

increment in 

I € 

. also semivertex ancle of delta 

J 2 „2 2 fl 2 2 

x+'V/x -p y -p z 

I 

plan form 

| ■ c 

2C 

- complex conical ray parameter 

l+€ 2 

I 

parameter defining r (defined differently in 
equations (l), (2), and (6)) 

m K 

constant 

J| n 

Mach angle (= sin~^ 


jj V 

pz/x 

1 - 

gas density 


1 • ' 

* 
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0 

T 

f 

Subscripts: 

A 

A,b 

A, c 
B 

B, b 
B,c 
i 

L.E. 

x 

7 

z 

1 

1, 2/3, 4 
I 
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Py/x except where designated otherwise 

dummy variable used in limiting process 

velocity potential; also. amplitude of elliptic 
integrals 

along OA (upper surface) 

along OA (upper surface), contribution of 
flow 7(b) 

along OA (upper surface), contribution of 
flow 7(c) 

along OB (right-band side) 

along OB (right-hand side), contribution of 
flow 7(b) 

along OB (right-hand side), contribution of 
flow 7(c) 

dummy index 

along leading edge 

dx 

_a_ 

dy 

dz 

r tinning variable 
distinguishing subscripts 
in region I 


II 


• in region II 
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APPENDIX ~B 


DELTA WINS WITH SYMMETRIC LINEAR TWIST 
(SUBSONIC LEADING EDGES) 

The loading on a delta wing with symmetric linear s panwise 
variation of angle of attack is easily obtained when the leading 
edges are supersonic. The result is given in equations. (16) and 
(17), upon multiplication by -pV. The difficulty is increased 
considerably when the leading edges are subsonic. The subsequent 
derivation applies to this subsonic case. 

The procedure employs the development of Busemann’s complex- 
variable method for conical flows. presented in reference 10. The 
storting point is an assumed function for the distribution of u 
with two undetermined parameters a and b. The procedure of 
reference 10 can then be applied to obtain a value of bv/bx in 
terms of a and b, and then again to obtain bv/by in terms of 
a and b. Insertion of the condition for linear twist (w = -ply I) 
requires that bv/bx = 0 and bv/by = ^P. The simultaneous solu- 
tion of these two equations for a and b completes the solution 
for u. 

The proper form for the assumed function for u was arrived 
at partly from the requirements of reference 10 and partly from a 
consideration of the asymptotic form assumed by equation (17) when 
the leading edges are sonic (m = 1). A limiting process for the 
sonic case yields 

u = 2 £X_ 1+jl. ( B1 ) 

3rtp 2 Vl“0 2 

Comparison of this result with the known u-velocity for a pitching 
delta wing for m = 1 and for m ^1 suggests the generalization 

u _ 2px a(m)+c£ b(m) _ (b2) 

3rtp 2 \jv?-o 2 

for the present wing for m£ 1, Other considerations, omitted 
for brevity, support the choice. Equation (B2) will reduce to (Bl) 
when m=lifa(l)=b(l)=l. 
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The relevant part of the conical-flow theory of reference 10 
is briefly as follows: It is proved that any function of the com- 
plex parameter 


r =- 2 £_ 


where 


fly+ipz 

x+<v/x 2 -3 2 y 2 -3 2 z 2 


gives a value of u that is conical and consistent with the 
Prandtl-Glauert equation (12). This value of u - call it U(£) - 
is complex in general , and both the real and the imaginary parts 
separately satisfy equation (12). Furthermore, the value of w 
associated with u = Be U may be obtained from the relation 


w(0 = - Re i3 




For z = 0 the parameter £ reduces to 3y/x = o. Thus by 
the foregoing principle a suitable known or assumed conical surface 
distribution of u,f(£) may be analytically continued off the sur- 
face as a spatial distribution by replacing a by £ and taking 
the real part. Either the surface or spatial value of w associ- 
ated with this distribution of u may then be obtained from 
equation (B4). 

In the present application, the requirement that u be coni- 
cal (function of a alone) is not met by equation (B2); the veloc- 


ity gradients u 


and u are, however, conical: 

u = 2 P am 2 - ( 2a+m 2 b ) o 2 

3 « 3 ^ (m 2 -tf *)^/ 2 


_ 2p (a+2m 2 b )o-bo^ 
3 " S (m 2 -o 2 ) 3/2 . 


JJHS 


NACA TN 2285 ' 

The variables u x , v x , w x may in this case replace u, v., v, 
respectively, in the method, of reference 10. Alternatively u y , 
v v may replace u, v, w. These substitutions are easily Jus- 

v w 

tified by differentiating the Prandtl-Glauert equation (12) and 
employing the irrotationality conditions (v x = u y , . . . ) . 

The analytic continuation of equations (B5) and (B6) is thus 
obtained by the substitution of £ for O : 


»*«) . He u x . Be 

3 « P 2 ( m 2 . ? 2) 3 / 2 


(B7) 


u~(£) ss Re U v = Re 

^ 7 o.o 3/2 


2p (a+2m^b) £-b (P 

(m 2 -S 2 ) : 

Also, the following equations correspond to equation (B4): 


(B8) 


w Y (0 = -Be ip 


Vy(C) = 


-Re ip 



(B9) 


(BIO) 


u y 


are seen to be wholly real 
Equations (B9) and (BIO) then 


In the plane £ = o , U x and 

in the ranges -m <0 < 0, 0 < O < m. 
show that w x and w y must be constant in these ranges . A dis- 
continuity in Vy may, however, occur in crossing (an infinitesi- 
mal distance above) the origin because the odd powers of £ in 
equation (B8) become imaginary there. This behavior is consistent 
with the desired boundary conditions already discussed: 


(Bll) 


H 

II 

O 

-m < 

<J < m 

H 

-m < 

a <0 


29 


= -P 


0 < 0 <m 


<f*/ r^ ,. /lyfc^K >' . - • !CS*7* *&-•.*'&.!? j -’ <yf 


NACA TO 2235 


It is convenient) and according to the foregoing discussion it. 
will suffice/ then, to evaluate equation (B9) for w x at o = m- 


and equation (BIO) for Wy at o = m- and at O = -m+, A limit- 


ing process must he used, including an approach from slightly 
above the y-axis (z/o) because of the singularity at a = 

(The procedure is mathematically equivalent to taking the "finite 
part" as done by Hadamard (reference 13).) Equation (B9) thus 
becomes 


w Y (m-) . Lim l Lim -Re ip 


0 \u — ►O 


m-T+iu 


6 ? 


(B12) 


and similarly for equation (BIO). 


The integrations are. most easily evaluated with the aid of 
the following elliptic-function substitutions, together with the 
tables of integrals in references 14 and 15: 


‘£1. 


k sn t 


1,2 2 . . 

y£ -m = k cn t 


C = dn t 


d£ = -k sn t cn t dt 


where 


t = F(0,k) 


0 = sin’ 1 

V 1-m 


(B13) 


k = a / 1-m 


The evaluation of equations (B9) and (BIO) by use of the limiting 
process of equation (B12) yields 


/ \ ~Bp 

w x ( m_ ' = 3np 


|(a+m 2 b)K(k) - (a+b)E(k)| (B14) 

1-m 2 L J 


w_(m-) 


-2p a+2m b n 


3« m' 


3 2 


(B15a) 
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and changing the upper limit to -m+T+iu 

*„(-») = + 2 

Jr 7 O 


nT 


(B15h) 


The left-hand aides of these three equations are now evalu- 
ated "by means of equations (Bll). - Equations (B15a) and (B15b) lead 
to the same equation, and there result but two simultaneous equa- 
tions for a and b: 


0 

3*0 i_ m 2 


J^a+m^b)K - (a+b)sj 


_ _ -2p a+2m z b Jt 

rv 


3Jt 


nr 


(B16) 


The solution of equations (E16) is , 


a(m) 


b = b(m) b 


3m' 


’(E-nnK) 


(l-2m 2 )E-Hn 2 K 

3m 5 (K-E) 

(l-2m 2 )E+m 2 K 




(B17) 


The values of a(m) and b(m) just obtained are to be sub- 
stituted in equation (B2). This completes the solution for the 
distribution of u- velocity (proportional to pressure) for the delta 
wing with linear symmetric twist and subsonic leading edges. 

A slight modification of the foregoing procedure yields the 
known pressure distribution for a pitching delta wing in a very 
simple manner. The boundary conditions (equation (Bll)) are merely 
replaced by w x = -p/p and Wy = 0. With these values, equa- 
tions (B14) and (B15) lead at once to the correct values of a 
and b for a wing pitching with angular velocity q = p/(J. These 
are the values a^(m) and b-]_(m) , respectively, obtained from ref- 
erence 12, expressions for which are given in appendix A. 

A natural question arises; namely, whether the procedure of 
this appendix could be used to solve the main problem of the rotat- 
ing cruciform with subsonic leading edges because the form of the 





mmm . 
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pressure-distribution function might be guessed from the supersonic- 
edge case in the limit as m — *1. If this method were capable of 
successful application, the exact solution rather than an approxi- 
mate one would be expected. Unfortunately, there are certain dif- 
ficulties. First, the pressure distribution is known for the 
other limiting case m — ►O, and the functional form is quite dif- 
ferent. Therefore the functional form for intermediate values 
of in cannot be guessed with confidence. Second, in the contem- 
plated application the analytic continuation of the pressure dis- 
tribution would require that U( £) possess an imaginary part, the 
evaluation of which may offer considerable difficulty. This can 
be seen as follows: Associate the function U^(£) with the hori- 

zontal panels and U 2 ( £) with the vertical panels. Call the cor- 
responding induced velocities at the horizontal panels w^(x,y) 
and w 2 (x,y), respectively. Then the boundary condition may be 
stated as Wi(x,y) + w 2 (x,y) = -py. Thus, because of w 2 , the 
gradients dwj_/dx and dw-j/dy will undoubtedly fail to be con- 
stants; this can occur only if U-jJO possesses a suitable imag- 
inary part for -m <_ £ = O S. m. (See equations (B9) and (BIO) and 
the discussion thereafter.) 
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TABLE I - VARIATION OF NA0/s 2 p WITH r/s FOR 
LOW-ASPECT-RATIO MULTIPLANAR DELTA-WING 
CONFIGURATIONS OF 3 AND 4 PANELS 


r/s 

0.3 

0.45 

0.6 

0.7 

0.75 

0.8 

0.85 

0.9 

0.95 

1.00 

Four panels 

NA0/s 2 p 

(cruciform) 

0.710 

1.176 

1.540 

1.672 

1.686 

1.657 

1.566 

1.384 

1.066 

0 

Three 

panels 

NA0/s 2 p 

0.671 

1.067 

1.298 

1.395 

1.385 

1.335 

1.235 

1.115 

0.851 

0 


For three panels the value of K at r/s = 0.85 given as 0.227 in 
reference 3 is clearly in error. Comparison with the approximate 
value 0.278 in reference 16 suggests that 0.272 was probably 
intended; the value NAJ2$/s 2 p = 1.235 above corresponds to K = 0.272. 
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Figure 5. - Boundary conditions for quadrant of rolling cruciform delta 
with subsonic leading edges. Section x = constant. 



Flow 6(a) Flow 6(b) 

(Fig.' 5, modified) 


Flow 6(c) 



Figure 6. - Superposition of simpler flows to obtain flow specified by modification 

of figure 5. 
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FlKure 12. - Contributions of the first two Iterations to -PC, for cruel 

form with subsonic leading edges. 







Variation of daroping-in-roll parameter with m for cruciform delta wing. Collected 
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